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Abstract 

We study the map which sends vectors of polynomials into their 
Wronski determinants. This defines a projection map of a Grassmann 
variety which we call a Wronski map. Our main result is computation 
of degrees of the real Wronski maps. Connections with real algebraic 
geometry and control theory are described. 



1 Introduction 



We study the map W which sends vectors of polynomials (/i,...,/p) to 
their Wronski determinants: 



W{h,...,f,) 



fi 
f'l 



fl 



p-i 



fp-i 

J V 



(1) 



Besides an intrinsic interest, this map is related to several questions of alge- 
braic geometry, combinatorics and control theory as we describe below. 

The following properties of the Wronski determinant are well-known and 
easy to prove: 

1. W{fi, . . . , /p) = if and only if /i, . . . , fp are linearly dependent. 

2. Multiplication of (/i, . . . , /p) by a constant matrix A of size p x p results 
in multiplication of W{fi . . . , fp) by detA. 

*Both authors are supported by NSF. 
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These properties suggest that our map W should be considered as a map 
from a Grassmannian to a projective space. We recall the relevant definitions. 

Let F be one of the fields R (real numbers) or C (complex numbers). 
For positive integers m and p we denote by Gf = G-p{m,m + p) the Grass- 
mannian, that is the set of all linear subspaces of dimension m in F'""'"^. 
Such subspaces can be described as row spaces of m x (m + p) matrices 
K of maximal rank. Two such matrices Ki and K2 define the same ele- 
ment of Gf if Ki = UK2, where U G GL{m,F). It is easy to see that 
GF{fn,m + p) is an algebraic manifold over F of dimension mp. We have 
GF(l,m+p) = FP™+''-\ the projective space over F of dimension m+p — 1. 

We may identify polynomials of degree at most m+p—1 in the domain 
of the map in (^Q) with vectors in F"'"'"'' using coefficients as coordinates, 
and similarly polynomials in the range of W with vectors in F™^"*"^ . Then, 
in view of the properties 1 and 2 of the Wronski determinant, equation 
will define a map Gf(p, "m + p) pp™-P. Alternatively, we can also identify 
polynomials of degree at most m + p—l with linear forms on F™"*"^ . Then p 
linearly independent forms define a subspace of dimension m in F'""'"^, and 
we obtain a map 

0: GF(m,m + j9) FP'^P, (2) 

which will be called a Wronski map. To understand the nature of this map, 
we use the Pliicker embedding of the Grassmannian. 

The Pliicker coordinates of a point in G-p{Tn.,m + p) represented by a 
matrix K are the full size minors of K . This defines an embedding of 

GF(m,m + p) to FP^, = ^ ^ - 1 . We usually identify Gf with 

its image under this embedding, which is called a Grassmann variety. It is a 
smooth algebraic variety in FP^. 

Let S C FP^ be a projective subspace disjoint from Gf , and dimF S = 
N — dimGF — 1. We consider the central projection tts : FP^\S pp™P^ 
and its restriction to Gf , 

05 = Ks\g. : Gf ^ FP"^f. (3) 

Then (ps is a finite regular map of projective varieties. When F = C this map 
has a degree, which can be defined in this the number of preimages 

of a generic point and is independent of S . This degree was computed by 
Schubert in 1886 (see [121111^] for modern treatment). 
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Theorem A When F = C, the degree of (ps is 



_ l!2!...(p-l)!(mp)! 
^""'^^ m\{m + l)\...{m + p-l)\- 



Projective duahty imphes that d{m,p) = d{p,m). Here are some values 
of d{m,p) 





m = 


2 


3 


4 


5 


6 


7 


8 


p = 


2 


2 


5 


14 


42 


132 


429 


1430 


p = 


3 




42 


462 


6006 


87516 


1385670 


23371634 


p = 


4 






24024 


1662804 


140229804 







p = 5 701149020 
In particular, 

1 / 2m \ 

d(m, 2) = , the m-th Catalan number. 

^ ' ^ m+1 \ m J ' 

The numbers d{m,p) have the following combinatorial interpretation: they 
count the Standard Young Tableaux (SYT) of rectangular shape p x m. 

To see that the Wronski map is a projection we choose a center Sq in 
the following way. Consider the p x [m + p) matrix of polynomials 



Eiz) 



I F{z) \ 
F'iz) 



(5) 



where F{z) = (2;"*+^ 2:^+p 2 . . . ^ 1). For a fixed z, the row space of 
this matrix represents the osculating [p — l)-subspace to the rational normal 
curve F : FP^ pp'ri+p-i ^^le point f[z). The space Poly^^ of all non- 
zero polynomials g G F[2;] of degree at most mp, up to proportionality, will 
be identified with FP'"^ (coefficients of polynomials serving as homogeneous 
coordinates) . 

We claim that the Wrosnki map (0) : Gp pp'^P can be defined by 
the formula 

K ^ <j>{K) = det ^ G Poly7, (6) 
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where is a matrix of size m x (m+p) representing a point in the Grassman- 
nian Gy ■ First of all it is clear that © indeed defines a map FP™?' ; 

changing K to UK, U G GL{m, F), will result in multiplication of the poly- 
nomial (f){K) by det U . Furthermore, this map when expressed in terms 
of Pliicker coordinates, coincides with the restriction to Gp of a projection 
of the form tts as in (j21), with some center which we call 5*0 . We do not 
need the explicit equations of 5*0, but they can be obtained by expanding 
the determinant in (jH)) with respect to the last m rows, and collecting the 
terms with equal powers of z. 

Now we verify that polynomial (piK) in (jH)) is a Wronski determinant. To 
see this, it is enough to consider the "big cell" X of the Grassmannian Gy , 
which is represented by the matrices K whose rightmost minor is different 
from zero. We can normalize K to make the rightmost m x m submatrix 
the unit matrix. If the remaining (leftmost) p columns of K are {kij), 1 < 
i < m, 1 < j < p, then 

^{K)=W{h,K,...,U,K), 

where 

fl,K{z) = Z"'+P-' - ki.iZ"'-^ - . . . - k^^u 
f2,K{z) = Z^'+P-^ - ki^2Z"'-^ - . . . - k^,2, 

fp,K{z) = — ki pZ^"" ^ — . . . — k^p. 

Coefficients of these polynomials correspond to p linear forms that define 
the row space of the matrix K = [{kij),I]. This proves our claim that (jH)) 
coincides with the Wronski map. 

For p = 2, this interpretation of the Wronski map as a projection is due 
to L. Goldberg |Hj. Her notation for Catalan numbers is different from our 
present notation. 

In this paper we study the real map 0, that is we set F = R. One 
motivation of this study is the following conjecture due to B. and M. Shapiro: 
If w & Polj^ is a polynomial all of whose roots are real, then the full 
preimage (^"^{w) of this polynomial consists of real points. In 0| we proved 
this conjecture in the first non-trivial case mm{m,p} = 2. On the other 
hand, when m and p are both even, there are polynomials w G Poly^^ 
which do not have real preimages under the Wronski map. So it was natural 
to ask the question, whether for some m and p one can give a lower estimate 
for the number of real preimages. To asnwer this question, we compute in 
this paper the topological degree of the real Wronski maps. 
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Notice the following important property of 0: it sends the big cell X 
of the Grassmannian into the big cell Y of the projective space consisting 
of those polynomials whose degree is exactly mp. Moreover, it sends the 
complement G-p\X into FP™^\F. When F = R these cells X and Y can 
be identified with R™^, in particular they are orientable, and the restriction 
of to X is a smooth map 

: X ^ F, 0((9X) C BY. (8) 

To define the degree of such map (see, for example ^5), we fix some orien- 
tations on X and Y . Then choose a regular value ?/ G F of 0, which exists 
by Sard's theorem, and define 

deg0 = ± ^ sgn det0'(x), (9) 

using local coordinates in X consistent with the chosen orientation of X, 
and any local coordinate at y . The degree deg / changes sign if one changes 
one of the orientations of X or F ; it is independent of the choice of local 
coordinates within the class defined by the chosen orientation of X, and of 
the regular value y . In Section 3 we will discuss a more general definition of 
degree which does not use special properties of the Wronski map and applies 
to all equidimensional projections of real Grassmann varieties. 

To state the main result of this paper, we need a definition. Consider 
the sequences a = (aj) of length mp whose entries are elements of the set 
{l,...,p}, and each element occurs exactly m times. Suppose that the 
following additional condition is satisfied: for every n G [l,mp] and every 
pair i < k from {1, . . . , p} , 

#{j G [l,n] : aj = i} > #{j G [l,n] : a, = k}. (10) 

Such sequences are called ballot sequences or lattice permutations ^^l- For 
given m and p, the set of all ballot sequences is denoted by Sm,p. There is 
a natural correspondence between T,m,p and the set of the standard Young 
tableaux of rectangular shape pxm [T7| Proposition 7.10.3]: we fill the shape 
with integers from 1 to mp putting one integer in each cell; if aj = i we put 
the integer j to the leftmost unoccupied place in the row i. (As usual, the 
row number increases downwards). 

Frobenius and MacMahon independently found that the cardinality of 
is d{m,p), the same number as in (jD), see for example, jini Sect. Ill, 
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Ch. V, 103] or T7", Proposition 7.21.6]. Of course, the coincidence of these 
numbers is not accidental [7| IT7j. 

Let o G Sm,p, 0" = A pair {pj.Ok) is called an inversion if j < k 

and (jj > CTfc. In terms of the SYT, and inversion occurs each time when for 
a pair of integers the greater integer of the pair stands in higher row than 
the smaller one. The total number of inversions in a is denoted by inva. 
Now we define 



I{m,p) 



E (- 



(11) 



It is clear that I{m,p) = I{p,m) because a pair of entries in a SYT is 
an inversion if and only if the same pair in the transposed SYT is not an 
inversion. This permits us to restrict to the case 

m>p>2 (12) 

in the computation of the numbers I{m,p). Recently, D. White |19^ proved 
that I{m,p) = iff m + p is even. For odd m + p satisfying ()12j). he found 
that I{m,p) coincides with the number of shifted standard Young tableaux 
(SSYT) of shape 

m+p—1 m+p—3 m — p + 3 m—p+1 



An explicit formula for the number of SSYT (see, for example, [11^ Proposi- 
tion 10.4]) gives I{m,p) = 

1!2! • ■ ■ (p - l)!(m - l)!(m - 2)! ■ ■ ■ (m - p + l)\{mp/2)\ 
(^_p + 2)!(m-p + 4)!---(m + p-2)! (^^^)! (^^)! ■ ■ ■ (^^^)!' 

when m + p is odd. SSYT appear in Schur's theory of projective represen- 
tations of symmetric groups, see, for example, p]. Here are some values of 
J(m,p) : 

m=34567 8 9 10 11 12 



p = 


2 


1 





2 





5 





14 





42 





p = 


3 





2 





12 





110 





1274 





17136 


p = 


4 







12 





286 





12376 





759696 




p = 


5 









286 





33592 





8320480 
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When p = 2 and m is odd, J(m, 2) = c/((m — l)/2,2), a Catalan number. 
The main resuh of this paper is 

Theorem 1 The degree of the real Wronski map (0) is ±I{m,p). 

Corollary 2 If m + p is odd then the real Wronski map (0) is surjective; a 
generic point y G RP™^ has at least I{m,p) real preimages. □ 

It follows from Theorem A, that for all y G Polj^^ , 

ca.Yd (j)^-^ (y) n G-R^m, m + p) < d{m,p). (13) 

This estimate is best possible for every m and p, see or the remark at 
the end of Section 2. 

If both p and m are even, then Theorem ^ gives deg (jR(m,m + p) = 0, 
and in fact in this case the preimage in Corollary |21 may be empty, as examples 
in [S] show. 

The lower bound in Corollary |21 with p = 2 is best possible, as the 
following example given in |1] shows: 

Example 3 For p = 2 and every odd m, there exist regular values y G 
RP^™ such that the cardinality of (j)~^{y) is I{m, 2) = d{{m - l)/2, 2). □ 

To each p-vector (/i, . . . , fp) of linearly independent polynomials one can 
associate a rational curve / = (/i : . . . : /p) in FP^~^ , whose image is not 
contained in any hyperplane. The following equivalence relation on the set 
of rational curves corresponds to the equivalence relation on the p-vectors of 
polynomials: 

f ^ g if f = i o g, where i is an automorphism of FP^~^. (14) 

If (/i, . . . , fp) is a coprime p-vector of polynomials, then the roots of 
W{fi, . . . , fp) coincide with finite inflection points of the curve /. Notice 
that Gr C Gc can be represented by p-vectors of real polynomials, and to 
each such p-vector corresponds a real curve /. When p = 2, / = /2//1 is a 
rational function. If the pair (/i,/2) is coprime, roots of W^(/i,/2) are the 
finite critical points of / . Thus our Theorem ^ has the following 

Corollary 4 Let X be a set of mp points in general position in C , sym- 
metric with respect to R . Then the number k of equivalence classes of real 



7 



rational curves in RP^~ of degree m + p—1 whose sets of inflection points 
coincide with X satisfies k > I{m,p) . In particular, for p = 2, this number 
k satisfies 

<k < d{m,2), if m is even, and (15) 
d{{m -l)/2, 2) <k <d{m,2), if m is odd. (16) 

Examples in 4 show that for every m, the lower estimates in (fT3j) and ()16p 
are best possible. So when m is odd, the Wronski map : G-R{m, m + 2) ^ 
RP^™ is surjective, while for even m it is not. □ 

We prove Theorem Q in Section 2. In Section 3 we discuss the definition 
of degree for arbitrary projections of real Grassmann varieties and interpret 
Theorem ^ in terms of control theory. 

For the case p = 2, the results of this paper were obtained in j^, with 
a different method based on We thank S. Fomin, Ch. Krattenthaler, F. 
Sottile and R. Stanley for helpful suggestions. 

2 Computation of degree 

In this section we prove Theorem ^ We fix integers m,p > 2. Consider 
vectors of integers k = {ki, . . . ,kp) satisfying 

< ki < k2 < . . . < kp < m + p, 

and vectors of real polynomials q = (gi, . . . , gp) of the form 

qi{z) = + ai,„_iz'™"i + . . . + ai,fc,2*=S 

q2iz) = ^-+1 + a2,„z"^ + . . . + a2,fc,z^^ (17) 

qp{z) = z^+P-' + ap,rr.+p~2z'^^P'' + ... + ap,k^z'^^. 

Suppose that all coefficients Ojj, ki<j<m + i — 2, 1 < i < p, are 
positive, all roots of the Wronskian = W{qi, . . . ,qp) belong to the semi- 
open interval (—1,0] C R, and those roots on the open interval (—1,0) are 
simple. The set of all such polynomial vectors q will be denoted by 6(k) . 
The greatest common factor of {gi, . . . , g^} is z'^'^ . 

It is easy to see that 6(k) parametrizes a subset of the big cell of the 
Grassmannian Gp{,(m, m + p): the representation of a point of G-R{m, m + p) 
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by a vector from 6(k) is unique. Setting ki = i — 1, 1 < i < p, we obtain an 
open subset 6(0, 1 . . . , p — 1) C Gn^m, m + p) . We define 

k = ki + {k2-l) + ... + {kp-p + l)>0. 

It is easy to see tliat k is tlie multiplicity of the root of PVq at for q G &(k) . 
Using coefficients of q as coordinates, we can identify 6(k) with a subset of 
j^mp-fc ^ introduce an orientation by ordering these coefficients: 

(18) 

In this sequence, coefficients of qj precede coefficients of qk for j < k, and 
coefficients of one polynomial qj are ordered according to their second sub- 
script. It is useful to place these coefficients into a Young diagram Y with p 
rows, such that coefficients of are in the i-th row, their second subscript 
decreasing left to right. For q G &(k), we denote the negative roots of the 
Wronskian W = Wq by 

In addition to these, there is a root of multiplicity k ai 0. We denote by Aq 
the Jacobi matrix of the map 6(k) — > Polj^~^, q ^ VFq, using coordinates 
(fT5|) in the domain and x = {xi,X2, ... ,Xmp-k) in the range, where —Xj 
are the negative roots of as in ()19|) . So the i-th row of this matrix Aq 
corresponds to Xi, and the j-th column to the j-th term of the sequence 
(fTH|) . When = 0, so that k = (0, 1, . . . ,p — 1) , and 6(k) is an open subset 
of G-£i{m,m + p) , we have Aq = 0'(q), the derivative of the Wronski map 
with respect to the chosen coordinates. 

For example, b{m — l,m + l,m + 2...,m + p — 1) consists of vectors 

q,iz) = z^ + ai,„_iZ™-\ q2iz) = z^^\ . . . , qp{z) = z^^~''\ (20) 

the Wronskian is 

W{z) = ip- 2)!^™^'-! Hp - l)z + p\a,,m-i) , 

and its only negative root is — p!ai^m-i/(p ~ !)• So 

det Aq = ]9!/(]9 — 1) > for q G &(m — 1, m + 1, . . . , m + p — 1). (21) 

This example will be later used as a base of induction. 
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We denote by E the set of all increasing homeomorphisms e : R>o — >■ 
R'>0) ^{t) < t for t > 0. Let n > 1 be an integer, and e & E. A thorn 
T{n,e) in R" is defined as 

{x = (xi, . . . , x„) e R!^o : < 1 < j < ^ - 1, x„ < e(l)}. (22) 

Notice that this definition depends on the ordering of coordinates in R". 
We always assume that this ordering corresponds to the increasing order of 
subscripts. 

Lemma 5 Intersection of any finite set of thorns in R" is a thorn in R" . 

Proof. Take the minimum of their defining functions e. □ 

Lemma 6 Let T = T{n,e) be a thorn in R"' = {(xi, . . . and U its 

neighborhood in H"''^^ = {{xo,Xi, . . . Then U'^ = 11(1 R>o^ contains 

a thorn T{n + 1, ei) . 

Proof. There exists a continuous function 6 : T R>o, such that f/+ 
contains the set {{xo,:k) : x G T, < Xq < 5(x)}. Let 6o{t) be the minimum 
of 6 on the compact subset {x G T(n, e/2) : Xi > t} of T. Then there exists 
eo E E with the property cq < Sq. If we define ei = min{e/2,eo}, then 
T(n + l,ei) C f/+. □ 

Lemma 7 Let T = T{n + 1, e) be a thorn m R"+^ , and h : T ^ R>^\ 
(xo,x) (?/o(xo, x), y(xo, x)) , a continuous map with the properties: for 
every x such that (xo,x) G T for some Xq > 0, the function xq ^ ?/o(xo,x) 
is increasing, and lim^^^o y(xo, x) = x. Then the image h{T) contains a 
thorn. 

Proof. We consider the region D G R""*"^ consisting of T, its reflection T' 
in the hyperplane Xq = and the interior with respect to this hyperplane of 
the common boundary of T and T' . The map h extends to T' by symmetry: 

Xo,x) = — /i(xo,x), (xo,x) G T, and then to the whole D by continuity. 
It is easy to see that the image of the extended map contains a neighborhood 
U of the intersection of D with the hyperplane xq = 0. This intersection is 
a thorn Ti in R" = {(xo,x) G R"+^ : Xq = 0}. Applying Lemma IHl to this 
thorn Ti , we conclude that contains a thorn. □ 

Given an increasing homeomorphism e E E , we define w{k,e) C Poly^^ 
as the set of all real monic polynomials of degree mp with mp — k negative 
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roots as in (|19p. these roots satisfying (j22|) with n = mp — k, and a root of 
multiphcity A; at 0. Thus w{k,e) is parametrized by a thorn T{mp — k,e) 
in BJ^P-'^. 

Starting with b{m — l,m + l,...,m + p— 1), we will generate subsets of 
6(k) by performing the following operations F*, 1 < i < p, whenever they 
are defined. Suppose that for some i G {1, . . . ,p} and some multiindex k, 
the following condition is satisfied: 

i > 1 and ki > ki_i + 1, or i = 1 and ki > 0. (23) 

Notice that for given k, this condition is satisfied with some i G {1, . . . ,p} 
iS k > 0. If (j^^ holds, we define a family of operators : 6(k) — * 6(k — ej) , 
where ej is the i-th standard basis vector in R^, by 

q^F^(q) = (q + a^'=-ie,), (24) 

where a > is a small parameter, whose range may depend on q. Thus 
an operation leaves all polynomials in q, except qi, unchanged. The 
following Proposition shows, among other things, that are well defined if 
the range of a is appropriately restricted. 

Proposition 8 Suppose that for some e E E , and k and i satisfying fj23|) . 
a set U C &(k) is given, such that the map q i— > W^^ : U — > w{k^ e) is 
surjective, and 

det Aq ^ for q G U. (25) 

Then there exist e* & E and a set U* C 6(k*), where k* = k — e^, with the 
following properties. Every q* G U* has the form F^{q) where F^ is defined 
m (j211), q G and a> 0; 

the map q* i— >• W^* : U* — » w{k — 1, e*) is surjective, (26) 
and det Aq. 7^ for q* G f/* . Moreover, 

sgn det Aq. = (-l)^(k,*) ggj^ (27) 

for every q* G U* and every q G t/, where x(k, «) «s t/ie number of terms in 
the sequence (|18p whose first subscript is less than i . In other words, x(k, i) 
is the total number of cells in the rows 1 to i — 1 in the Young diagram Y 
described after (fTSj) . 
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Proof. Let us fix q G U, and put W = Wq. As W E w{k,e), we have 
ord W = k, where ord denotes the multiphcity of a root at . Let cz^ be 
the term of the smallest degree in W{z). Then c > 0, because all roots of 
W are non-positive. In fact, 

c = l[{kj-k)l[a,,k,>0. (28) 

j>l 3 

We define W* = W^*, where q* = -^'(q). Then ordW* = k — 1 and the 
term of the smallest degree in W*{z) is c*2;^~^, where 

c* = al[ik*-k*)l[a,,k, >0. (29) 

We conclude that when a is small enough (depending on q), the Wronskian 
W* has one simple root in a neighborhood of each negative root of W , and 
in addition, one simple negative root close to zero, and a root of multiplicity 
— 1 at . To make this more precise, we denote the negative roots of W 
and W* by 

-Xn< ■■■ < -xi and - yn < . . . < -yi < -yo, (30) 

where n = 2m — k, and yj = yj{a) . We have 

yj{0) = Xj, for I < j < n, and yo{0) = 0. (31) 

Furthermore, if a is small enough (depending on q) 

a \~>- yo{a) is increasing and continuous. (32) 

The set w{k,e) is parametrized by a thorn T = T{n,e), where x = 
{xi, . . . , Xn) , and n = mp — k. There exists a continuous function 6o : T —>■ 
R>o, such that 

q*G6(k*), for aG(0,5o(x)), x G T. (33) 

Now we are going to compare det Aq with det Aq* . For this purpose we 
investigate the asymptotic behavior of the 'new root' yo{ct) of W* , as a — * 0. 
Comparison of the terms of the lowest degrees in W{z) and W*{z) , (j28|) and 
(EHl) show that 

- yo(a) = -c*/c + o(a) = -c^^a/ai^k, + o{a), a ^ 0, (34) 
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where Ck > depends only on the muhiindex k. 

The Jacobi matrix A* = Aq* is obtained from the Jacobi matrix A = Aq 
by adding the top row, corresponding to yo , and a column, corresponding to 
tti^ki-i = CL- The position of the added column is 

1 + X(k,0, 

where x(k, i) is the number of terms of the sequence (fTH| whose first subscript 
is less than i. 

According to (jHlI), the intersection of the added row with the added col- 
umn contains the only essential element of this row: 

dyo/da = c^/ai^k, + o{l), a 0. 

The rest of the elements of the first row of A* are o(l) as a — >■ 0. Expanding 
A* with respect to its first row, we obtain 

det A* = (-l)^(M (ck/ai,fcj det A + o(l), a ^ 0. 

Now it follows from our assumption that for sufficiently small a , A* 7^ . 
Moreover, ()27|) holds, if a is sufficiently small. More precisely, for every 
q G t/ there exists 61 (q) > such that for < a < 5i (q) we have det A* 7^ , 
and (|77j). Taking 6 = m.m{6o,6i} , where 60 was defined in (jH^ . we obtain 
the set 

U* = {q* = : X G T, a G (0, 5(x))} C b{k*), (35) 

and this set U* satisfies (P7|) . Here qx G t/ is some preimage of G 
w{k, e) = T under the map q Wq. Such preimage exists by assumption of 
Proposition IHl that the map q 1-^ Wq, U w{k,e) is surjective. It remains 
to achieve (|26|) by modifying the thorn T. This we do in two steps. First 
we apply Lemma IHl to the half-neighborhood (|35|) of T, with xq = a, to 
obtain a thorn Ti(n+ l,ei) in R""*"^. Then we apply Lemma [7| to the map 
/i : Ti — >• R>o^, defined by yj = yj{xo,x.), where yj are as in (jHUj) . and 
Xq = a. This map h satisfies all conditions of Lemma [7| in view of ()31|) and 
(jSai- This proves n 

Conclusion of the proof of Theorem Q| We begin with a brief outline 
of our argument. For each k and i satisfying (|^. equation defines 
an operator depending on parameter a: : 6(k) 6(k — ej). Starting 
from a subset of b{m — 1, m + 1, . . . , m + p — 1) , we will consecutively apply 
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operators in all possible sequences allowed by (j23j) . In the end we will 
obtain a set of polynomial p-vectors in 6(0, 1, . . . ,p — 1) , which will contain 
the full preimage of a point under the Wronski map. Equations fl27|) will 
permit to control the sign of the Jacobian determinant of the Wronski map 
at all points of this preimage. 

Now we give the details. Consider the set of all finite (non-empty) se- 
quences a = (cTj), where cXj G {!,... ,p}, j G N, satisfying (fTUI) . For every 
such sequence we define 

k(cr) = {ki, . . . , kp), where ki = m + i — 1 — #{j : aj = i}, 



Let S = S(m,p) be the set of all sequences a satisfying (fTUI) and k{a) > 0. 
Notice that for k = k((T), condition (P^jl holds with some i G {1, . . . ,p} if 
and only if k{a) > 0. 

To each sequence a G S we put into correspondence an open set 11^ C 
6(k(cr)) in the following way. For cr = (1), we set 

Uii) = {qG 6(m- l,m + 1) : G u;(l,eo)}, 

where eo(a;) = x. Then f/(i) consists of the polynomial vectors of the form 
with and ai,„_i G (0, {p-l)/p\). 

Applying operations F* to f/(i) means that we use Proposition |H1 with 
U = [/(I) , and k = (m — 1, m + 1, . . . ,p— 1) . We obtain from this Proposition 
the sets U* , which we call f/(i,i) . In fact, Proposition |H1 can we applied 
in this situation only with i = 1 or i = 2. Then we apply operations 

, j G {1, . . . ,p} to , whenever permitted by and so on. 

In general, suppose 11^ is already constructed. If k(o") and i satisfy (j^ . 
we apply operation to Ua ■ This means that we use Proposition |S1 with 
U = Ua, k = k(cr) and this i. The resulting U* is called t/((T,i) C 6(k(cr)— Cj). 

Every sequence a G S encodes an admissible sequence of applications of 
operations F*. If a = (cTj), then aj = i indicates that F* was applied on 
the j-th step. Conditions (fTUI) and k{a) > imply with some i, so that 
an operation F* is applicable. Every operation decreases k{a) by 1, so the 
procedure stops when k{a) = 0. 



and 



p 




i=l 
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Proposition ISl implies that for each cr G S with k{a) > 0, there exists 
& E such that 

W ■U,^w{k{a),e,) (36) 

is surjective and unramified. 

Observe that we can always replace e* in Proposition |H1 by a smaller func- 
tion from the set E . We use this observation to arrange that the coefficient, 
added to polynomials in q on each step, is strictly smaller than all coeffi- 
cients added on the previous steps. This implies that for each q G f/^, all 
coefficients are strictly ordered, and the sequence a can be recovered from 
this order. More precisely, let k = k{a) , and ci > C2 > . . . > C2m-k > 
be the ordering of the sequence of coefficients of qi . . . ,qp. Then aj = z if 
Cj = 0,1,1 with some /. In other words, enumerating the cells of the Young 
diagram Y defined after (jl8|) in the order of decrease of their entries gives a 
standard Young tableau. The sequence a can be recovered from this tableau 
in a unique way. 

We recall that the number of inversions inva was defined in the intro- 
duction, just before the equation (fTT|) . We claim that for every cr G E, 

sgndetAq = ^(a)(-l)^°"" if q G f/,, (37) 

where /i(cr) = ±1 depends only on the length of a. Indeed, by (j27j), on each 
step the sign of det A is multiplied by (—1)'^, where x = x(k(cr),i) is the 
number of terms of a which are less than i. This proves (jHTj) . 
Now we consider the subset 

Sm,p = {o- G Em : k{a) = 0}. 

It consists of ballot sequences, as defined in the introduction. The set S^.p 
corresponds to rectangular standard Young tableaux of the shape p x m. 
The number of such tableaux is d{m, p) (see, for example, f7| Proposition 
7.21.6]). Sequences a G Sm,p generate d{m,p) open sets Ua- C 6(0, 1, . . . ,p — 
1) with the property that the maps ()36p are surjective and unramified. Using 
Lemma IHl we restrict these maps so that they have a common range w{0, e) 
with some e G -E. 

As all maps (j36p are surjective, every point from this common range 
has at least one preimage under the Wronski map in each U^,a & ^m,p- AH 
these d{m,p) preimages are different as elements of 6(0, 1, . . . ,p—l) , because 
the sequence a can be recovered from the ordered sequence of coefficients 
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of q G [/o-. Furthermore, all these d{m,p) polynomial vectors represent 
different points in the Grassmannian G-[i{m,m + p) , because to each point 
in 6(0, 1, ... ,p — 1) corresponds only one point of G-R^m, m + p) . Thus we 
found d{m,p) different preimages of a point under the Wronski map. On the 
other hand, the complex Wronski map has degree d{m,p) by Theorem A, so 
we found all preimages of the real or complex Wronski map. Equation (j37|) 
gives the signs of Jacobian determinants at these points, so the degree of the 
Wronski map is given by □ 

Remark. In the process of this proof, we constructed a point in RP™^ 
which has d{m,p) distinct real preimages under the Wronski map. This 
proves the fact earlier established by Sottile [ini, that the upper estimate 
d{m,p) given by (fT^ . is best possible for every m and p. 



3 Additional comments 

1. Let us show how to define topological degree (an unsigned integer) for 
arbitrary projections of real Grassmann varieties as in (jS)) with F = R. 

Let f : X ^ Y he a. smooth map of compact, connected real manifolds 
of equal dimensions. If X is orientable, the degree deg / can be defined by 
formula Q. If ^ is orientable but Y is not then deg / = 0. 

Now we suppose that both X and Y are non-orientable and consider 
canonical orientable 2-to-l coverings X ^ X and Y ^Y , which are called 
the spaces of orientations of X and Y jH 10.2]. The set X consists of pairs 
{x, O) where x G X and O is one of the two orientations of the tangent space 
Tx . There is a unique structure of smooth manifold on X which makes the 
map X ^ X, (x, O) (— s> X a covering, and O depends continuously on x. 
The group of the covering X — > X is {±1} . Notice that the spaces X and Y 
have canonical orientations. A map / : X — >■ F is called orientable if there 
exists a lifting f : X ^ Y , which commutes with the action of {±1}, d 
(10.2.5)]. A different but equivalent definition of an orientable map is given 
in [13 §5]. 

For an orientable map, we define deg / := ± deg /. Under our assumption 
that X and Y are connected, this degree is defined up to sign, which depends 
on the choice of the lifting. Though Y is connected, Y may consist of one or 
two components, but the degree is independent of the choice of the regular 
value y eY . 
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Suppose that there exists a regular value y ^ Y and an affine chart 
f/ C X, so that f~^{y) C U. Then we can compute the sum (jU)) using 
coordinates in U . The orientability of / ensures that this sum is independent 
of the choice of the chart U and coincides with deg / . 

To apply this construction to projections of Grassmann varieties 

we recall that RP™^ is orientable iff mp is odd, and G-£i{m,m + p) is ori- 
entable iff m + p is even (see, for example, [HI Ch. 3 §2]). So in the case that 
m + p is even, the degree of (ps is defined in the usual sense, as in Q. 

To deal with the case when m + p is odd, we first identify the space of 
orientations of a Grassmannian G-Ei{m,n), with odd n. Consider the "upper 
Grassmannian" G^{m, n) , which consists of all oriented m-subspaces in R". 
It can be also described as the set of all m x n matrices K of maximal rank, 
modulo the following equivalence relation: K' ^ K if K' = UK, where 
detU > 0. We have the natural 2-to-l covering G^{ 
which assigns to the class of K in G^{m,n) the class of the same K in 
GR(m, n) . We also have n) = (RP"~^)+ , a sphere of dimension n — 1 . 

Every upper Grassmannian is orientable and has canonical orientation. 
To see this, we consider the tangent space = Tx{G^{m,n)) which is 
the product of m copies of a subspace y = R"""*, complementary to x. 
Orientation of x induces a unique orientation of y, such that x®y = R" has 
the standard orientation. This defines an orientation on each tangent space 
Tx which varies continuously with x. So we have a canonical orientation of 
(^^(m, n) . 

We claim that for odd n, the coverings 

G^{m,n) G-R_{m,n) and GR,(m, n) — » (^^(m, n) are isomorphic. (38) 

Indeed, for x G (^^(m, n), orientation of x C R" defines an orientation 
of Tx{G-£i{m,n)) , as explained above. One can easily show that (in the 
case of odd n) changing the orientation of x changes the orientation of 
Tx{GR{m,n)) . This proves (jHHj) . 

We recall that a projection map tts '■ RP^\S' RP'^ can be described 
in homogeneous coordinates as 

y = Ax, (39) 
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where A is a (k + l) x {N + 1) matrix of maximal rank, and x, y are column 
vectors of homogeneous coordinates. The null space of A represents the 
center of projection S = S{A) (where the map is undefined). Two matrices 
define the same projection if they are proportional. A change of homogeneous 
coordinates in RP^ or in the target space RP'^ results in multiplication of 
A by a non-degenerate matrix from the right or left, respectively. 

Proposition 9 Let n he an odd integer, G-R,{ra^ n) C RP^ a Grassmann 
variety, and : GR,(m, n) p^pm(n-m) ^ central projection. Then is 
orientable. 

Proof. The Pliicker embedding PI : GR,(m, n) RP^ lifts to G^{m,n) - 
(RP^)+, which is defined by the same rule as PI. Using (|38|). we identify 
Gii{m,n) with G^{m,n), and obtain the lifting 

PI : GR(m, n) (RP^)+ (40) 

of PI. Suppose now that a projection vr^ is defined by fl39|) where A is an 
{m,{n — m) + 1) x [N + 1) matrix. Then the same equation (j39|) defines the 
lifting 

7r+ : (RP^)+\5+ ^ (Rp™("-"^))+ ^ pp"^("-™)^ (4^) 

where is the preimage of S under the covering (RP^)+ — > RP^, and 
the last isomorphism holds because m{n — m) is even. Composition of the 
maps pUj) and is the desired lifting of (ps , which is evidently compatible 
with the action of {±1}. The existence of such a lifting proves that (ps is 
orientable. □ 

Thus equidimensional projections of real Grassmann varieties always have 
well-defined degrees. It is clear that when the center of projection varies 
continuously, the degree does not change until the center S intersects the 
Grassmann variety. These exceptional centers form a subvariety Z of codi- 
mension 1 in the Grassmannian Gb.{N + 1,N — mp) of all centers. So the 
degree is constant on every component of G'r(A^ + 1, — mp)\Z , in partic- 
ular, all projections (ps whose centers S belong to the same component of 
G-R,{N + 1, — mp)\Z as the center 5*0 of the Wronski map have the same 
degree ±I{m,p). 

2. Now we restate our results in terms of control theory by static output 
feedback. Suppose that a triple of real matrices S = {A, B, G) of sizes 
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n X n X m and p x n is given. This triple E defines a linear system 

X = Ax ~\~ BzL^ //ir)\ 
y = Cx. ^^"^^ 

Here x, u and y are functions of time (a real variable) taking their values 
in R", R*" and R^, respectively. The values of these functions at a point 
t G R are interpreted as the state, input and output of our system at the 
moment t. 

Behavior of the system (j42j) is completely determined by its transfer func- 
tion z ^ C{zl — A)~^B, which is a function of a complex variable z with 
values in the set of p x m matrices. One wishes to control a given system 
by arranging a feedback, which means sending the output to the input 
via an m X p matrix K , called a gain matrix: 

u = Ky. (43) 

Elimination of u and y from (j42j) . (j43p gives the closed loop system 

x = {A- BKC)x, 

whose transfer function has poles at the zeros of the polynomial 

ijK (z) = det {zI~A- BKC) . (44) 

The map K i— »• ipx ^ PoIYr is called the pole placement map, and the 
problem of pole assignment is: given a system S, and a set {zi, . . . , Zn} , 
symmetric with respect to R, to find a real gain matrix K , such that the 
zeros of ipK c^'^e {zi,...,Zn}- Thus for a fixed system E, arbitrary pole 
assignment is possible iff the pole placement map is surjective. 

When n > mp, X. Wang ^H] proved that for generic S, the pole place- 
ment map is surjective. Here we consider the case n = mp. We also as- 
sume that n is the smallest possible size of a matrix A in the representation 
C{zl — A)^^B of the transfer function. Systems with this property are called 
"controllable and observable", and and they form an open dense subset of 
the set of all systems with fixed (m, n,p) . To understand the structure of the 
pole placement map, we use a coprime factorization of the open loop transfer 
function of a generic system S (see, for example, j21 Assertion 22.6]): 

(45) 
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C{zl - A)-^B = D{z)-^N{z), det D{z) = det{zl - A), 

where D and are polynomial matrix-functions of sizes p x p and p x m, 
respectively. The polynomial matrix [D{z), N{z)] has the following prop- 
erties: its full size minors have no common zeros, and exactly one of these 
minors, det D{z), has degree n while all other minors have strictly smaller 
degree. Every p x {m + p) polynomial matrix with these properties is related 
to a linear system of the form (j42j) via equations (j45j) . 

Using the factorization and the identity det{I — PQ) = det{I — QP) , 
which is true for all rectangular matrices of appropriate dimensions, we write 

ijK{z) = det{zl -A-BKC) = det{zI -A)det{I -{zl -A)-^BKC) 
= det{zl - A)det{I -C{zl - A)-^BK) 
= det D{z) det(/ - D{z)-^N{z)K) = det{D{z) - N{z)K). 



This can be rewritten as 



D{z) N{z) 
K I 



Poly7. (46) 



In the last determinant, the first p rows depend only on the given system, 
and the last m rows on the gain matrix. Permitting arbitrary m x {m + p) 
matrices K of maximal rank as the last m rows of the determinant in ()46p 
we extend the pole placement map to 

0s : GR(m, m + p)^ RP"^^ ct>^{k) = [^kI (47) 

where [.] means the class of proportionality of a polynomial, which is identi- 
fied with a point in RP*"^ , using the coefficients of a polynomial as homoge- 
neous coordinates. The map (j47p is defined if for every matrix K of rank m 
in the last m rows of ()46|) the determinant in pUj) does not vanish identically. 
Systems S with this property are called non- degenerate, and they form an 
open dense subset in the set of all systems with given {m,p) and n = mp. 
Applying Laplace's expansion along the first p rows to the determinant in 
(|46p ■ we conclude that the map 0s , when expressed in Pliicker coordinates, is 
nothing but a projection of the Grassmann variety G-£i{m,m+p) into RP*"^ 
from some center depending on S . This interpretation of the pole placement 
map as a projection comes from Now we notice that all projections 

arising from linear systems as in pUj) have the property that they send the 
big cell X of G-[i{m, m + p) represented by matrices K of the form [K, I] 
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into the big cell Y of Poly^^ consisting of polynomials of exact degree mp. 
Furthermore, GR,(m,m + p)\X corresponds to Poly^^\y inder such projec- 
tions. Our arguments in the first part of this section imply that the real pole 
placement maps of a non-degenerate system has a well-defined degree^. As 
the center of projection S = varies continuously, this degree remains 

constant as long as S does not intersect the Grassmann variety. Degenerate 
systems are precisely those for which S intersects the Grassmann variety. 

Comparing (jH)) with ()47p we conclude that the Wronski map is a pole 
placement map for some special linear system. So our Corollary El implies 

Corollary 10 For every m and p such that m + p is odd there is an open 
set U of linear systems with m inputs, p outputs and state of dimension mp, 
such that for systems in U the real pole placement map is surjective. Fur- 
thermore, the pole placement problem for systems in U has at least I{m,p) 
real solutions for any generic set of mp poles symmetric with respect to the 
real line. □ 
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